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Abstract 

We analyze the classical stability of string cosmologies driven by 
the dynamics of orientifold planes. These models are related to time- 
dependent orbifolds, and resolve the orbifold singularities which are 
otherwise problematic by introducing orientifold planes. In particu- 
lar, we show that the instability discussed by Horowitz and Polchinski 
for pure orbifold models is resolved by the presence of the orientifolds. 
Moreover, we discuss the issue of stability of the cosmological Cauchy 
horizon, and we shoiu that it is stable to small perturbations due to 
in-falling matter. 
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1 Introduction 



Over the past year there has been a renewed interest in the study of time- 
dependent string backgrounds [1-33] . The main objective of this programme 
is to learn about the origin of the (possibly apparent) cosmological sin- 
gularity using the powerful techniques of string theory. The logic behind 
this investigation was to start with an exact conformal field theory, e.g. 
strings in flat space, and to consider an orbifold of the theory with a time- 
dependent quotient space "^i^. These simple orbifold constructions were, 
however, quite generally shown to be unstable by Horowitz and Polchinski 
[TB] ^. These instabilit ies are associated to the formation of large black holes 
whenever a particle is coupled to the geometry and are non-perturbative in 
the string coupling. For this reason a perturbative string resolution of the 
time-dependent orbifold singularities seems inappropriate. 

In IH] we investigated the orbifold of flat space-time by a boost and trans- 
lation transformation. This orbifold suffers from the Horowitz-Polchinski 
instability. However, an aditional essential ingredient of the model was in- 
troduced in 1^: the time-like orbifold singularity is an orientifold, i.e. a 
hrane with negative tension. One of the results that we shall show is that 
the orientifold provides precisely the non-perturbative string resolution of 
the singularity, modifying the gravitational interaction in its proximity. In 
fact, after briefly revising the Horowitz-Polchinski argument and the orien- 
tifold cosmology construction, we shall show in Section 4, with a tractable 
example in dimension 3, that due to the presence of the orientifolds large 
black holes are not formed. 

In any pre big-bang scenario, there is generically a possible classic in- 
stability associated with the propagation of matter through the bounce, 
which is analogous to the propagation of matter between the outer and in- 
ner horizons of charged black-holes [361 137j . Let us consider a perturbation 
at some finite time in the far past, when the universe is contracting. In the 
case of orientifold cosmology there is a future cosmological horizon where 
the contracting phase ends (see also [SHI jBHl EDI)) therefore one needs to 
worry about such perturbations diverging at the horizon, thus creating a 
big crunch space-time. We shall show in Section 5 that such perturbations 
remain finite (of course, they can grow and create small black holes, but this 
does not change the causal cosmological structure of the geometry). Fur- 

^The exception is the nuU-brane |35J with extra non-compact directions |14|. which is 
an exact and regular bounce in string theory free of singularities. Notice, however, that 
this space is not a cosmology since it is supersymmetric and depends on a null, not a 
time-like, direction. 
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thermore, we shall see that the fluctuations that destabilize the geometry 
are precisely those that already destabilize the Minkowski vacuum. 

2 The Horowitz— Polchinski Problem 

In ^H] , Horowitz and Polchinski argue that a large class of time dependent 
orbifolds are unstable to small perturbations, due to large backreactions of 
the geometry. Their results, which overlap with the string theory scattering 
computations of [3 E]) do not rely on string theory arguments, and are 
obtained simply within the framework of classical general relativity. 

The argument of Horowitz and Polchinski is quite simple and runs as 
follows. They analyze the orbifold of D-dimensional Minkowski space 
by the action of a discrete isometry e'^ where k is a Killing vector, and 
consider, in the quotient orbifold space WP /e'^ , a small perturbation. This 
perturbation naturally corresponds to an infinite sequence of perturbations 
in the covering space M^, related to each other by the action of e'^. They 
consider, quite generally, a light ray with world-line Q. C M^, together with 
all of its images 

0„ = e"^r^ (n G Z) 

for any integer n. The basic claim of '16" is that, for a wide range of time 
dependent orbifolds, the backreaction of gravity to the presence of this in- 
finite sequence of light rays 0^ will inevitably produce large black holes. 
As we already mentioned, this statement is a claim in pure classical gravity 
theory and to prove it Horowitz and Polchinski consider the scattering of 
two light rays and They compute, given k, the impact parameter 
hn and the center of mass energy of the scattering, and show that, for 
generic time dependent orbifolds, one has that G8n ^ b^~^ for n sufficiently 
large. The two image particles will be well within the Schwarzchild radius 
corresponding to the center of mass energy, and will form a black hole. 

It is quite clear that, if the above argument is correct, it must hold also 
for a continuous distribution of light rays 

e^^n {t G M) (1) 

since the problematic interactions occur between rays separated by At ^ 1, 
and therefore should be completely insensitive to the distribution of the 
essentially parallel light rays when the parameter t varies by At ~ 1. For 
continuous values of t, equation defines a surface S in M^, which we call 
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the Horowitz-Polchinski (HP) surface. Again, the backreaction of gravity 
to the matter distribution on S should form black holes by the previous 
argument. On the other hand, the continuous distribution is easier to treat 
analytically since it is invariant under the action of k. Therefore, the full 
solution with the gravitational backreaction will also have Killing 
vector. This fact simplifies explicit computations and therefore, from now 
on, we will only consider the continuous case. 

In the rest of the paper we will be interested exclusively in a specific 
choice of orbifold, which describes, when correctly interpreted in M-theory, 
the dynamics of an orientifold/anti-orientifold pair. We postpone the discus- 
sion of the correct M-theory embedding and the crucial distinction between 
the orbifold and the orientifold interpretation to section |31 and we describe 
now the orbifold along the lines of jl6j . 

We write = x E^-^ ^^^^ 

we parametrize with coordinates 

X^, Y, with line element 

-dX^dX- + . (2) 

Points in M'^ are identified along the orbits of the Killing vector field k which 
corresponds to boost in the X^ plane and translation in the orthogonal 
direction Y 

K = X+d+ - X-d- + dy . (3) 

For this special choice of k one has bn ^ n and £n ~ e"" and the bound 
G£n S> b!^~^ is always satisfied for n large ^Hl- The directions along E^~^ 
then play no role (except to determine the strength of the gravitational 
interaction) and we will omit any reference to them, working effectively in 
three dimensions. 

To study the geometry of the HP surface S, we parametrize the light ray 
Q as 

A" + sV , 

where V"" is a null vector and s is an affine parameter. If we consider 
the generic case^ when V^V^ ^ then, for some value of t, the image 
ray e*'^0 will be directed along the null vector^ (1, 1, 1). This specific light 
ray on the surface S will intersect the plane X'^ + X~ = at a point 
with X~^ — X~ = 2a. By shifting the variables s,t,Y so that this point 

^When V^V~ — the images e'Ti correspond to parallel light rays, and do not form 
black holes. 

It could also be directed along (1,1,-1). This case is just the mirror under X"^ 
and we do not discuss it. 
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Figure 1: The curves given by equation ©, which represent the surfaces Sa 
for the values of a indicated in the figure. The surface 5*1 is clearly singled 
out, and is the unique null surface. 

corresponds tos = t = Y = Owe see that the surface S is given by the 
following parametrization 

X± = {s±a)e^\ (4) 
Y = (s + t) . 

The only free parameter which labels different surfaces is the constant a, 
and we will add a label Sa to distinguish them. A convenient way to visualize 
Sa is to move to coordinates 

x± = X±e^^, (5) 

y = Y, 

so that K = Then, since the surface 5"^ is invariant under k, it is given 
by a curve in the x^,x~ plane (see Figure ^) 

x^ = (s±a)e^^ (6) 

This curve can also be viewed as the intersection of the surface Sa with the 
plane y = in M^. 
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Let us now describe some basic properties of the surfaces Sa- The in- 
duced metric is given by 

2 (1 - a) dsdt + (s^ - + l) dt'^ 

so it is timehke for all values of a except for a = 1, when 5" is a null surface. 
For o 7^ 1 the surfaces are smooth for any value of the parameters s,t. For 
a = 1 the surface is, on the other hand, divided into two smooth parts for 
s > and s < which are joined along the singular line s = 0. Moreover, 
the forward directed light cones, which are tangent to the null surface Si, 
lie on one side of 5i for s > and jump discontinuously, along the singular 
line s = 0, to the other side of Si for s < 0. 

2.1 The Three— Dimensional Horowitz— Polchinski Problem 

As we described in the previous section, we want to understand, in general, 
the gravitational backreaction to a distribution of light -rays in x E^~^ 
distributed on the surface Sa parametrized in by ^ and fixed at a point 
along the spectator directions E^^^. For D > A gravity is non-trivial even 
in the absence of matter, and the problem is not soluble exactly. In D = 3, 
on the other hand, space is curved only in the presence of matter and the 
problem is tractable. 

Recall that the basic question is whether black holes are formed as a 
consequence of the strong gravitational interactions among the light rays. 
On the other hand, black holes in three dimensional flat space do not exist 
in the first place, and therefore it seems that the question is not well posed 
for D = 3. This is not correct. In fact, the sign of the formation of a 
black hole in three dimensions is not the presence of a singularity but the 
existence, in the geometry, of closed timelike curves. To see this, recall 
again that, in a two-particle scattering with center of mass energy £ and 
impact parameter b, a black hole forms when G£ > b^^^. For D = 3 
the threshold is then G£ > 1 and is independent of the impact parameter 
b. To understand the meaning of this threshold, recall that the geometry 
describing the propagation of a single light ray in 3D is flat everywhere 
except along a null line, where curvature is concentrated. The local geometry 
is completely characterized by the holonomy in 5*0(1,2) around the light 
ray, which is given by a specific null boost Ui. In the presence of a second 
light ray, we have a second holonomy null boost matrix U2, as shown in 
Figure 121 The combined holonomy U = U1U2, on the other hand, depends 
on the center of mass energy of the two particles. For small center of mass 
energies U is a rotation. On the other hand, when the center of mass energy 
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\J=U,U, 
CTC 



Figure 2: Scattering of light rays in three dimensions. The holonomies Ui 
and U2 are null boosts, whereas {7 is a boost for center of mass energies 
greater then a specific threshold. In this case, the geometry has CTC's, as 
indicated. 

exceeds the threshold value [ISl El ESI IS] j U becomes a boost and closed 
timelike curves are formed which circle both particles, as shown in Figure 
121 Finally, as shown in [IJ, one may consider the same dynamics in AdS^ 
space, where black holes exist I47J . One then finds that, at exactly the same 
threshold, the scattering process creates, in the future, a BTZ black hole. 

We now see that the three dimensional problem consists in analyzing 
the gravitational backreaction to light rays located on the surface Q), and 
investigating the existence, in the full geometry, of closed timelike curves. 
As we already mentioned, since the surface Sa preserves the Killing vector k, 
the full geometry with backreaction will also have the same Killing vector. 
We will show two things: 

1. The geometry will have closed timelike curves. It is important to stress 
that we are talking about the uncompactified geometry, and not of the 
one obtained after identifying by e'^. 

2. If we excise from the spacetime manifold the region where k is timelike, 
then the resulting geometry is free from closed timelike curves. As 
discussed in [S], the excision is necessary to correctly embed these 
geometries in M-theory. This procedure introduces a boundary in the 
geometry which corresponds, in string/ M~theory, to the presence of 
orientifold planes. We will recall these points more at length in the 
next section. 
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3 Orientifold versus Orbifold Cosmology 



As we just mentioned, and as discussed in detail in j^, when we embed 
the orbifold spaces in string/M-theory we must excise the regions with 
k2 < 0. The space is naturally divided into two regions, X and y, 
where the vector k is spacelike (or null) and timelike, respectively. Since 
= 1 + X~^X~ , X and y correspond to the regions X^X~ > —1 and 
X^X^ < —1. We can now consider the quotient manifold 

m3 = MVe-^ = XQUyQ, 

where Xq = X je^ and yq = y /e^ . The action of k, is free and therefore the 
quotient manifolds are smooth. On the other hand, the full quotient space 
Mq, although regular and without boundary, has closed timelike curves 
which always pass in region 3^q. The space Xq, on the other hand, does not 
have closed timelike curves, but does have a boundary. 

In the usual time dependent orbifold model, one considers the following 
compactification of M and II A theory 

X 

/si K \ 

M| X ^ ? 

The top line is the M-theory vacuum. Going to the left, we are compacti- 
fying along a circle of the 8-torus, and we obtain the Type IIA orbifold 
with constant dilaton. Going to the right, we are compactifying along the 
Killing vector k. The bottom right corner should be related, as usual, by 
a TST duality transformation to the Type IIA orbifold on the bottom left. 
On the other hand, the IIA supergravity background on the right is quite 
peculiar, since it has a complex dilaton field. This is because the vector k 
is timelike in yq. 

In jni EI we have analyzed more carefully the bottom right corner of 
the diagram and have proposed a related but different compactification of 
M-theory, which uses crucially non-perturbative objects of the underlying 
string theory. In fact, it was noted in that, if one excises from the M- 
theory vacuum the region yq and dimensionally reduces the space Xq x 
along K, one obtains a warped IIA supergravity solution with non-trivial 
dilaton and RR 1-form of the form Z xT^, where Z is a two-dimensional 
space with boundary. The boundary of Z is singular and corresponds to 
the supergravity fields of an orientifold plane, which acts as a boundary of 
spacetime, and has well defined boundary conditions for the string fields. 
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The geometry is best understood after T-dualizing along the 8~torus and 
describes then the interaction of a 08-08 pair, a system which has also 
been studied in conformal field theory in [481 149j . We now have a consistent 
compactification scheme given by 

Xq X <^ Z xT^ 

where we are adding tildes on the spaces to recall that one has to impose 
specific boundary conditions on the string fields at the boundaries of these 
spaces. 



4 Reducing to Two Dimensions 

We now move back to the analysis of the three-dimensional Horowitz- 
Polchinski problem. As we already mentioned, the full geometry after back- 
reaction will still have a Killing vector. Therefore, upon dimensional reduc- 
tion, the problem is effectively two-dimensional. 

Let then the general form of the metric in three dimensions be 

dsl = dsl + $2 {dy + Aadx'^f , (7) 

where ^ is the Killing direction. The three-dimensional Hilbert action 
reduces to 



The equation of motion for the field strength F is trivial and implies that 
the scalar field (^^-kF is constant. By rescaling the variables y, Aa and 
we can fix the constant to any desired value (provided it does not vanish) 
so that 

Then, the equations of motion for the dilaton and the metric can be derived 
from the action 



Jd'xVg(^R-^y 



(9) 



This action is a particular case of two-dimensional dilaton gravity (see |5Ul 
I51| for reviews and a complete list of references), and we will find it most 
convenient to work in this framework from now on. Only at the end we will 
connect back to the discussion in three dimensions. 
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We conclude this section by showing that Q is also relevant in the study 
of the interaction of 08-08-planes wrapped on a 8-torus, which again is 
naturally a two-dimensional problem. In general, the system is described 
by the following massive IIA background 



E'dsio = ^ dsi + ds^r) , (10) 

2E 
9s 



-5 ^ 2^ 



Then the Type IIA action 

becomes the simple two-dimensional dilaton-gravity action 



d^x^ <!>R 



which again reduces to © after rescaling of ^. Note that, in conformal field 
theory, we can construct only a single 08-08 pair with specific tensions and 
charges. This then implies (see .9 for more details) that 

Els = 2gs. (11) 



4.1 Two Dimensional Dilaton Gravity 

In this section we review some basic facts about 2D gravity which will be 
useful in the analysis of the action Q . 

Two-dimensional dilaton gravity is a natural way to define, in two di- 
mensions, theories with a gravitational sector. The basic fields are the mat- 
ter fields, as well as the metric gab and a dilaton field which should be 
considered together as the gravitational fields. The action has the general 
form 

S2D ((7,^') + 5m (5,^, Matter). 

Sm is the matter part of the action, whereas 6*20 is the generalization of the 
Einstein-Hilbert term and is given by 

S2D= I d^x^[^R-Vm , 
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where V {^) is a potential for the dilaton. The equations of motion are 
easily derived to be 



2VaVt,$ = gab (2n$ + V)- Tab , (12) 



dV 



where Tab and p are 



2 dSu 1 (^"Sm 



Moreover, the conservation of the stress-energy tensor Tab is modified by 
the dilaton current p to 

+ pVb^ = 0. 

The inherent simplicity of the dilaton gravity model lies in the following 
observations jHSEOlEl]- Define J ($) by 

J = J Vd<l> 

and consider the function 

C7=(V$)VJ($) (13) 

and the vector field 

= —e'^'^Vb^ . 

Then, whenever Tab = p = - i.e. for any vacuum solution to the equations 
of motion - the function C is constant and k is a Killing vector of the 
solution. The first fact follows immediately from the equations of motion 
which imply 

VaC = -TabV'<^ + Va<^ {Tbcd"') ■ (M) 

The second fact is proved most easily in conformal coordinates , with 
metric 

-dz+dz'e^ . 

Then we have that 
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and the non-trivial Killing equations become V+V+<I> = V_V_<I> = 0, 
which now follow from (|12() whenever Tab = 0. Finally note that these 
equations are equivalent to 

d-K~^ = d-^K~ = . 

Using these facts it is trivial to find all classical vacuum solutions. In 
fact, whenever < 0, the metric can be put locally into the form 

d 

ds^ = - f {x) dt^ + dx^ , K = a—, (15) 

where a is a constant. Then, one has the explicit solution of the equations 
of motion 

which depends only on the constant C. Similar equations hold in regions 
where > 0. 

Let us now analyze the geometry in the presence of matter. For our 
purposes, we are going to consider only matter lagrangians Sm which do 
not depend on the dilaton, and which are conformal. This implies that 

T+- = p = 0, 
9_r++ = 5+r = . 

The simplest example is clearly a conformally coupled scalar r/ with action 
Sm = — J (V??)^. The effect of this type of matter is best described by con- 
sidering a shock wave (42 1 illj , which is represented in conformal coordinates 
by a stress-energy tensor of the form 

T.^{z-)=e6{z- -z,) . (e>0) 

The positivity of e can be understood by looking at the conformally coupled 
scalar, for which r = 2 (V„r/)^ > 0. Recalling from (|14p that 

V_C = 2r__V+$e"^ = r__K-, (17) 

we conclude that the shock front interpolates, as we move along z~ , between 
the vacuum solution with C = Cq and the vacuum solution with C = Co + 
eK~ {zq) (see Figure El. As a consistency check note that, since in the 

vacuum r and k.~ are functions only of z~ , equation ()17() defines a jump 

in the function C which is independent of the position z^ along the shock 
wave. 
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Figure 3: Shock wave solution in two-dimensional dilaton gravity. 



4.2 The Case V = 2$-^ 

We now specialize to the action Q by taking 

_ 2 _ 1 

First consider the vacuum solutions. These are easy to obtain by considering 
the flat metric © written in the coordinates = X^e^^^ and y = \fEY ^ 
where > is a constant with units of energy, 

^ = Ve - x-d) + ^dY 

dy ^ ' ^ 

is a Killing direction and are the two-dimensional coordinates. If one 
rewrites the metric in the form ((7|, it is easy to check the normalization Q 
and to compute the two-dimensional metric and dilaton fleld. The explicit 
computation is done in [3] and the result corresponds to the cosmological 
solution shown in Figure lU where the timelike orientifold singularities are 
at a distance of order . For future use, we divide the diagram in various 
regions lin,out, ^^L,R and 111l,r, as shown in the Figure. 

Analytic results are simplest in Lorentzian polar coordinates. For exam- 
ple, in the Rindler wedges 11/, and II/j, where x^x^ < 0, one uses coordinates 
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Figure 4: The cosmological vacuum solution for V = 2<I>~^. Regions lin, 
lout are the contracting and expanding cosmologies, whereas regions IIl,r 
are the intermediate regions. Regions 111l,r are after the singularities, and 
correspond, when uplifted to three dimensions, to the regions where the 
Killing vector k becomes timelike. 



t, X defined by x = ±xe and obtains the static solution ^ 

„2 



dsl 



x 



' 1 - E^x^ 



+ dx^ 



(18) 



It is easy to check that the above solution satisfies (|15|) and pfij) with 

C = -E. 

As we already discussed in section |31 these C < solutions correspond to 
non BPS 08-08 geometries when uplifted to string theory with equation 
(|1U() . Notice that, given E {oi equivalently Qs by equation CJ)), one has no 
freedom in the solution, which is unique. Therefore, no fine tuning is required 
in the initial conditions for the metric and the dilaton in order to obtain a 
solution with a bounce and with past and future cosmological horizons. Recall 
also that the conformal field theory description of the 00 system has two 
distinct moduli, which are the string coupling and the distance between 
the orientifolds in string units. We therefore see that the backreaction of 
the geometry ties the two moduli with equation Ul\) and reduces them to a 
single free parameter. 
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We also notice that the C = solution is clearly singled out and corre- 
sponds to the supergravity solution of a single BPS OS-plane. The C > 
solutions do not have a clear interpretation in string theory, and it would 
be nice to understand this point further. 

Before considering more general solutions in the presence of matter, let 
us discuss the coordinate change from to conformal coordinates z^. For 
the sake of simplicity, and since this is all we will need, we write formulae 
for the case E = 1. We consider then the change of coordinates 

+ - 2$ ^ - 1 x+ 

= e ^ 7, — = — , (19) 

$ + 1 z- X- ^ ^ 

where 

$ = \/x+x" + 1. 
Then, it is simple to show that the metric becomes conformal 

, 9 dz+dz- _i\ 

dsl = -^TI- ( ] ■ (20) 



Z^Z' 



Finally let us add matter and consider shock wave solutions. Given 
the discussion in section (|4.1() . it is now almost trivial to find the correct 
geometries, which are given pictorially in Figure [51 In Figure before the 
shock wave, we have the vacuum solution with some given value of E. After 
the wave, one has again a vacuum solution, but with a different value E' . 
Recall that 

E' = E- eK- 

where e > and = 2V+<I> e^^ must be computed along the wave. As 
we move in Figure from point a to points b and c along the shock wave, 
in the direction of increasing z^ , the value of $ decreases to at c on the 
singularity. Therefore V+$ < and one has that 

E' > E. 

Moreover, in any vacuum solution with C = —E, the value of the dilaton on 
the horizons is l/\/^, as can be seen from equation (|18() at a; = 0. Therefore, 
since the value of ^ is continuous across the shock wave, we notice that the 
horizon to the left of the wave, where the dilaton has value ^ = l/\/^, 
must intersect the wave between the points b and c, as drawn. 

A particular case of the shock wave geometry is attained when the wave 
moves along one of the horizons, as shown in Figure 13). Since <I> is constant 
along the horizons, we have that V+<1> = and therefore that 

E' = E. 
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a 

(a) 



Figure 5: Shock wave solutions in the cosmological geometry. 

Let us breifly explain why the horizon at a constant value of z'^ shifts as 
one passes the shock wave. It is easy to see that the horizon in question is 
given by the curve = 0. In the vacuum, is a function of alone, but 
in the presence of matter one has that 

d^K+ = e-^T__ . 

Then, since is constant along the horizons (and therefore along the shock 

wave) and since r has a delta singularity, the function k+ just jumps by 

a finite constant across the wave, thus explaining the shift in the position of 
the horizon. 

Let us conclude by noticing that the instability of the null orbifold found 
in is again a shock wave which interpolates between two different vacuum 
solutions. The initial solution is a BPS vacuum (C = 0), which has a 
different spacetime structure from the other solutions with C 7^ which 
are formed after the shock wave. In our case this instability does not arise 
since we are interpolating between two non-BPS vacua with the same global 
structure. 
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Figure 6: The shock wave before backreaction, both in conformal coordinates 
and in the coordinates of equation H19() . 

4.3 The Shock Wave Solution in Three Dimensions 

The geometry of Figure [5^ corresponds, when uphfted to three dimensions, 
to the gravitational backreaction to matter distributed along a specific sur- 
face, which is the uplift of the shock wave curve. In order to better under- 
stand the geometry of this surface, consider first the shock wave without 
backreaction. From now on we take, for simplicity, E = 1 and E' = 1 + 6 
for some 6 > 0. The limiting geometry for 5 — > is given by Figure El The 
underlying 2D geometry is the vacuum E = 1 solution, and we have singled 
out a specific curve, which is given, in the conformal coordinates z^, by the 
two segments 

z- = -1, (-oo<z+<l) (21) 

z+ = 1 . (-1 < < oo) 

We now uplift this configuration to three dimensions. The vacuum solu- 
tion becomes, as we emphasized before, nothing but flat space. The curve 
becomes then a specific surface in M^, whose geometry is most easily ana- 
lyzed by passing to the coordinates x^ in (fT9|) . It is simple to see, using the 
explicit coordinate transformation ()19() . that the curve H21|) becomes (see 
Figure 

x^ = {s± 1) e^^ (22) 

for s G M. The first part of the curve (PT|) corresponds to (|^ for s < 0, 
whereas the second part corresponds to s > 0. Comparing with Q we 
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immediately see that this curve corresponds, when uphfted to three dimen- 
sions, to the Horowitz-Polchinski surface ^i. Recall from section |21 that, 
among the various surfaces Sa, the surface Si is special since it is the only 
null surface, and it is therefore not surprising that we are obtaining exactly 
Si from the two-dimensional shock wave. 

We therefore conclude that the geometry of Figure when uplifted 
to three dimensions, represents the gravitational backreaction to light rays 
distributed on the HP surface Si . 

4.4 Finding the Closed Timelike Curve 

We are now in a position to investigate the existence of closed timelike curves 
in the geometry of Figure IHK- Recall that, without shock wave, the quotient 
geometry has closed timelike curves that always pass through region III, and 
that are removed by excising this region. These curves are not closed in the 
covering space, which is flat. In the geometry with the shock wave, we shall 
see that there are closed timelike curves already in the covering space, which 
could signal an instability. However these curves always pass through region 
III, and again should be excised. 

In order to describe the closed timelike curves, we first need to prove 
a simple fact. Consider, in M^, two points A and B which are connected 
by a future-directed timelike curve from A to B. Let x^,0 and x^,yB 
be the coordinates of these two points in the coordinate system ^ and 
let us assume that both points are in Region II^j, i.e. the region with 
x+ > 0, < and x^x^ < 1. We parametrize the path from A to B as 
x^ (s) ,y (s), with < s < 1. First of all, it is easy to show that the full path 
must be entirely regions II/j or III_r, because otherwise the horizons in the 
geometry prevent the curve from returning to region IIjj without becoming 
spacelike. Therefore, since we are confined to the right Rindler wedge of the 
x^ plane, we may adopt polar coordinates = itxe^*. We want to show 
that, if tA > ts (the point A is after the point B in Rindler time), then 
any forward directed timelike curve from A to B must go into region III r 
(the region with x > 1). To prove this fact, recall first the metric on in 
the coordinates x, t, y 

dx^ - x^dt^ + dy"^ (l - x^) - Ix^dydt . 

Consider then the function t [s). It starts at t (0) = t^ and it is increasing 
for small values of s. This is because surfaces of constant t are spacelike 
in region 11r (since dx"^ + dy"^ (l — x^) > if x^ < 1) and therefore future 
directed timelike curves have always 4^ > in region II/j. Now, since by 
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assumption ts < t^, the function t{s) must achieve a maximum for some 
value of s. Then, using the same reasoning, the curve must be in region 
IIIr to be timehke, and this concludes the proof. We can actually always 
construct a future directed timelike curve by connecting the two points A 
and B with a straight line, provided that we choose i/b large enough, as can 
be seen by direct computation. 

We can now prove the existence of a closed timelike curve in the full 
geometry (jS^)- The regions to the left and to the right of the wave are 
vacuum solutions which correspond, when uplifted to three dimensions, to 
regions in flat Minkowski space. We can then think of the full geometry 
as follows. Let and be two copies of three-dimensional space which 
are going to describe the geometry to the left and to the right of the shock 
wave, and which we parametrize with coordinates 

= x^e^ Y = y (M^) 

respectively. The shock wave itself defines surfaces S and S in the two 
spaces and M^. The full geometry is then obtained as follows. Denote 
the region to the left of S in by M|. This corresponds to the part of 
the geometry to the left of the shock wave. Similarly, denote the region 
to the right of S in by M|j, which is the region to the right of the 
wave. The full geometry is then given by taking the two regions M^^ and 
M|j and gluing the two boundaries S and S" in a specific way. We do not 
need the details of the gluing, aside from the simple fact that y = y when 
connecting S and S. Looking at Figure [7J which represents the shock wave 
solution in conformal coordinates, choose then the two points A and B as 
drawn. They will correspond to points on the surface S with coordinates 
x^,yA and x^,yB and on S with coordinates x^,yA and x^,yB, where, as 
we already noticed, one has yA = UA and jjB = Ub- As in the discussion 
in the previous part of this section, we choose yA = and ye large and 
positive. Notice first that, with respect to the geometry on the right of 
the shock wave, both points A and B are in region 11^, whereas they are in 
regions lout and Ij„ relative to the geometry M.^ on the left of the wave. In 
equations, we have that 

^ 0, A ^ ^' A*^ A ^ ^' 

B ^ B ^ B^ B ^ ' 
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Figure 7: A closed timelike curve in the geometry induced by the surface 
Si. 

and 

> 0, 
S| < 0. 

This is the key feature which is required in order to have a closed timelike 
curve. Let us first concentrate on the right part of the wave. We may just 
consider the straight line from yl to i? in Mjj which is future directed for ys 
large, and which goes, when projected onto the plane as in Figure [3 in 
the region after the singularity. This is because the point A is after point B 
in the Rindler time of region IIr (x^/xg > x\/x^), and so we can apply 
the arguments of the first part of the section. Now we consider the region 
to the left of the shock wave. It is easy to check, by direct computation, 
that the straight segment from to ^ in is still timelike for large values 
of y_B = Ub, and it is clearly future directed since B is in Ij„ and A is in 
lout- This closes the loop, as drawn in Figure [71 and we have found a closed 
timelike curve. 

Note that, due to the simple lemma proved above, it must be that our 
closed curve always goes into region III. Recall though that the boundary 
of region III corresponds in string theory to the orientifold 08-plane, and 
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Figure 8: A closed timelike curve in the geometry induced by the surface Sa 
for a < 1. 

therefore region III should be excised, as described in section IHI The final 
M -theory geometry is free of closed timelike curves, and we have then shown 
that the presence of the orientifolds not only resolves the issues put forward 
in section\^ but also cures the instability due to the formation of large black 
holes. 

4.5 Discussion 

Let us conclude this section with some comments. 

First of all, we can try to consider the geometry induced by light rays 
on surfaces 5a for —1 < a < 1 (so that the surface does not go into region 
III). This corresponds, in 2D gravity, to matter which couples to the dilaton 
and to the conformal factor of the metric, and is therefore analytically more 
complex. Nonetheless, we can understand pictorially, in Figure |H1 that the 
physics is qualitatively unaltered. In particular one can find, as before, 
closed timelike curves which necessarily have to pass in region III. 

Secondly, it is clear that, in FigureEt) it is crucial that E' > E. This fact, 
recall, comes from positivity requirements on the stress-energy tensor Tab- If 
E' had been less then the graph would have had a different structure. In 
particular, the horizon of the left part of the geometry would have intersected 
the shock wave not between the points b and c, but between a and b. The 
causal structure of the space, when uplifted to three dimensions, is then 
quite different and it is easy to show, using arguments similar to those in 
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the previous section, that the 3D geometry has no closed timehke curves. 
Finally, let us comment on the general Horowitz-Polchinski problem for 

D > 3. In this case, gravity is non-trivial even in the absence of matter, 
and therefore an exact solution to the problem is probably out of reach. On 
the other hand, it seems unlikely that a correct guess on the final qualita- 
tive features of the solution can be obtained by looking at the interaction 
between two (or, for that matter, a finite number) of light -rays. This fact is 
already true if we just consider the linear reaction of the gravitational field 
to the matter distributed on the HP surfaces Sa- Then, very much like in 
electromagnctism, it is incorrect to guess the qualitative features of fields by 
looking at just a finite subset of the charges (matter in this case), whenever 
the charge distribution is infinite (this infinity is really not an approxima- 
tion in this case, since it comes from the infinite extent of the surface Sa 
due to the unwrapping of the quotient space). Therefore, to decide if the 
problem exists in higher dimensions, much more work is required, already 
in the linear regime of gravity, but most importantly in the full non-linear 
setting. Note that the only case in which the HP argument is fully correct 
is exactly in dimension D = 3, where the gravitational interaction is topo- 
logical and when, therefore, the interaction of an infinite number of charges 
can be consistently analyzed by breaking it down into finite subsets. This 
indeed is what we find in the previous sections, since closed timelike curves 
do appear. 

5 Stability of the Cosmological Cauchy Horizon 

Another related classical instability which can arise in orbifold and orien- 
tifold cosmologies is due to the backreaction of general matter fields as they 
propagate through the bounce: as the universe contracts, particles will ac- 
celerate and will create a large backreaction in the geometry. General insta- 
bilities of this type can only be studied in the linearized regime, as opposed 
to the shock wave geometries which could be analyzed fully, including grav- 
itational backreaction. In fact, in the previous section we concentrated on 
conformally coupled non-dilatonic matter propagating in the cosmological 
geometry, and this was possible due to the simple coupling to the two- 
dimensional gravitational fields $ and gab- We now want to consider more 
general matter fields, either non-conformal or coupled to the dilaton. In the 
following, we shall concentrate on a scalar propagating in the quotient space 
M^/e'^, which corresponds to a free massive scalar on the covering space M^. 
Our results will apply to the orientifold cosmology after imposing Dirichlet 
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or Neumann boundary conditions on the surface X^X^ = —E"^. These 
wave functions were studied in [9|, and we shall review briefly those results 
before considering the issue of stability. 

5.1 Particle States 

Let us start with a massive field ^ of mass m on the covering space M'^, 
which satisfies the Klein-Gordon equation 

U^{X) = rr?^{X) . (23) 

We demand that ^ be invariant under the orbifold action, so that 

^[X) = ^{e'^X) , (24) 

where k = 27rz (A J + RK) and 

iJ = X+d+ -X-d- , 
iK = dy . 

The energy scale E of section |1] is given by -E = AR^^. Since the quotient 
space has the symmetry generated by the Killing vector k, and since the 
operators J, K and □ commute, it is convenient to choose a basis of solutions 
to and (Plj) where the operators J and K are diagonal. We then choose 
our field "if so that 

J^{X) = p^{X) , 
K'^iX) = k'^iX) , 

by writing 

(X) = ^p(X+,X-)e*^^. 

The wave functions are functions of Lorentzian spin J = p which also 
satisfy the two-dimensional Klein-Gordon equation (49+9_ + a;^) '^p = 0, 
where u;^ = m? + k"^. They can then be solved in Lorentzian polar coordi- 
nates, as we shall review below, in terms of Bessel functions, which naturally 
have an integral representation as a sum of the standard plane waves in the 
covering space. In terms of the quantum numbers p, fc, the orbifold bound- 
ary condition (|24() becomes simply 

Ap+ Rk = neZ. (25) 

Generally, when R ^ 0, the spectrum of p E M is continuous for each re. On 
the other hand, for the pure boost orbifold {R = 0), the spectrum is discrete 
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since Ap = n. This difference is crucial. We shall see more in detail in the 
next section that, in the general case, it is possible to define particle states 
that do not destabilize the cosmological vacuum on the horizons by properly 
integrating over the various spins p, whereas in the pure boost case there 
is an infinite blue shift as particles approach the tip of the Minkowski cone 
j34j . The situation is analogous to particles propagating in the "null boost" 
and "null boost + translation" orbifolds studied in [14| . From now on we 
consider the general case R ^ 0, and we use as quantum numbers p and n, 
thus labeling the wave functions ^'p^. 

To study the behaviour of the wave functions describing the collapse 
of matter through the contracting region Ij„ we first move to light-cone 
coordinates = X^e^^^ and then to polar coordinates in the Milne wedge 
of the plane x^ = te^^^. This means that, using the quantization 
condition ()25|). the wave function ^'p,n is given by 



and where J±ip are the Bessel functions of imaginary order zizip. In terms 
of the quantum numbers p and n, the frequency to satisfies the mass shell 
condition 



5.2 Near Horizon Behaviour of the Wave Functions 

We are interested in the limit t ^ near the cosmological horizon, so it is 
convenient to recall the basic functional form of the Bessel functions 



where F±ip is an entire function on the complex plane, which has the ex- 
pansion 



Then, in terms of the two-dimensional light-cone coordinates x^, which are 
well defined throughout the whole geometry, we have 




(26) 




•x 
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These wave functions are well behaved everywhere except at the horizons 
= 0, where there is an infinite blue-shift of the frequency. To see this, 
consider the leading behaviour of the wave function "^^^ at both horizons 
(x"*" = or x~ = 0) 

' \ 2 J r{i + ip)' 

Near the horizon x~ = the wave function is well behaved and can be 
trivially continued through the horizon. Near x~^ = 0, on the other hand, 
the wave function has a singularity which can be problematic. In fact, close 
to the horizon, the derivative of the field d^"^^^ oc (x"*")*^ ^ diverges as 
x+ ^ 0, and this signals an infinity energy density, since the metric near 
the horizon has the regular form ds"^ ~ —dx^dx^. This fact was noted 
already in pilin|. 

A natural way to cure the problem is to consider wave functions which 
are given by linear superpositions of the above basic solutions with different 
values of p. The problem is then to understand if general perturbations in the 
far past t <C —E~^ will evolve into the future and create an infinite energy 
density on the horizon, thus destabilizing the geometry. This problem is well 
known in the physics of black holes where, generically, Cauchy horizons are 
unstable to small perturbations of the geometry j36U37j . We will show in the 
next section that the problem does not arise in this cosmological geometry. 
More precisely, we will show that perturbations which are localized in the 
X direction in the far past t <C —E~^ (and which are therefore necessarily 
a superposition of our basic solutions (|26() ) do not induce infinite energy 
densities on the horizons of the geometry, and can be continued smoothly 
into the other regions IIl,_r and lout of the geometry. 

This problem was also considered in [21] , where a different conclusion was 
reached. We will discuss, at the end of section the argument of and 
describe why it does not apply to physically relevant bounded perturbations. 

5.3 Analysis of Infalling Matter 

For simplicity we will set, from now on, £^ = 1 and m = and we will 
focus uniquely on the uncharged sector n = 0, which corresponds to pure 
dimensional reduction. The three-dimensional action for the massless scalar 
^ reduces, in two dimensions, to 

j d^xy/g^ {V^f . 
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The scalar field ^ is therefore confer mally coupled, but has a non trivial 
coupling to the dilaton, and satisfies the equation of motion 

D^- + Vln$ • = 0. (27) 

As in the previous section, we analyze this equation in the background 
described by the vacuum solution (fTH)) . which is given in regions lin,out by 

^2 

ds'^ = -dt^ + ^dx^ , (28) 

1 + 

$ = Vl + t2 ^ 

where = te^^ . We have seen, in the previous section, that the solutions to 
(|27j) are known exactly in terms of Bessel functions, since (|27|) is equivalent 
to the equation (Dpiat — <^^) ^' = 0, where J = —i^ is the boost operator 
in the plane and Dpiat is the Laplacian with respect to the flat metric 
—dt^ + t^dx'^ . A different and more general approach to the solution of l|27() . 
which is more in line with the literature on the stability of Cauchy horizons 
in black hole geometries, is to rewrite H27() in terms of the field 

A = ^V^. 



A simple computation shows that equation H27|) becomes 

1 



□A + 



4$2 



A = 0, 



We now use the fact that, in the vacuum solution (|28j) . equations (fT2j) and 
(fni) with C = -1 imply that □^> = -2^'^ and that (V^>)^ = -1 + ^>"^. 
Therefore the equation of motion becomes 

5- $2 

As customary, we write the above equation using conformal coordinates 
defined in region In the sequel, we follow closely the beautiful work of 
Chandrasekar and Hartle i37| . We first define, in terms of the coordinates 

in (|19|1 . the coordinates u, v given by 

z+ = -e-\ 
z- = -e-"". 
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Figure 9: The tortoise time coordinate s versus t, together with the asymp- 
totic behaviours t and — ln(— t) for t — > — oo and t ~ 0. 

The coordinates u, v are hghtcone coordinates in the s, x plane 

u = s + X , 
V = s — X , 

where s is the usual tortoise time coordinate. Recalling s is given in 
terms of t by the expression 



1 VTT^ + i 



.._VTT15^_,„^__ ,30) 
and has asymptotic behaviour for large values of |s| given by (see Figure IH]) 

s ~ -ln(-t) . (s » 1) 

The various coordinate systems are recalled in Figure 112b . The metric 
in region Ij„ is given, from equation (|20j) . by —dudvt^ (l + t^) , and we 
conclude that equation (|29j) reads 

5i2 -a^)^(^'^) = ^(«)^(^'^) ' 



where 

^2 (4-^2) 

4 (l + t2) 

Finally, we Fourier transform the x coordinate 



y = - \ . ■ (31) 
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Figure 10: The scattering potential V{s). 



Ais,x) = / dpe'P''A{s,p) 



to obtain the Schrodinger-like equation 

(^-^+p^^A{s,p) = V{s)A{s,p) . (32) 

Let us consider this differential equation in more detail. The potential 
y (s), shown in Figure ITUl has the following asymptotic behaviour 

V ~ e-^' , (s^oo) (33) 

V ~ , (s^-oo) 

and is therefore significant only in the region around s ~ 0. The scattering 
potential connects the two regions with s <C — 1 and s ^ 1 (or t <C — 1 
and t ~ 0) where the field A behaves essentially like a free scalar in two- 
dimensional flat space. The potential V then connects the past Minkowski 
region (t <C — 1 with metric —dt"^ + dx"^) to the future Milne wedge (t ~ 
with metric —dt'^ + t^dx^). Since V decays at infinity faster then one 
can follow the usual theory of one-dimensional scattering. In particular, 
following we consider the solutions F{s,p) and P{s,p) of which 
behave like pure exponentials e~*^'^ respectively in the future s — > oo and in 
the past s — > — oo 

F{s,p) ~ e-'P' i-tyP , (s^oo) (34) 
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Therefore, when restoring the x dependence e*^^ of the fuh solution, we 
have, for s — > oo, the plane waves 



and, for s — > — oo. 



ipv 



As explained in the previous section, the problem is exactly soluble in terms 
of Bessel functions^, and one has the following explicit form of the functions 
F and P 



F{s,p) = - r(i + ip) (i + t2)4 



p) 

where t = t{s) is defined by ((SOJ. The function 

Hi'J = — rV^ (e'^'^.p - J-^p) (35) 



is the Hankel function of the first kind, which has the correct asymptotic 
behaviour to ensure that P{s,p) has the desired properties. 

It is useful, again following closely 0^, to consider the general analytic 
behaviour of F and P as a function of the complex variable p, at fixed time 
s. First it is clear that we must have 



F{s,p) = F {s, -p) , 
P{s,p) = P{s,-p). 

Moreover, we recall from that, whenever the potential has asymptotic 
behaviour 

V ~ (s^oo) 

V ~ e^'^+' , (s ^ -oo) 



^Recall that (z) = z" x (entire function of z). In this paper we choose unconven- 
tionally to take the branch cut of the logarithm along the negative imaginary axis, so that 
2" and (2) also have a cut on Re 2 = 0, Im2 < 0. Moreover, the expressions like p^^ 
and will always be defined using the same prescription for the logarithm. 
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i n 



F(s,p) 



P(s,p) 



Figure 11: The analytic structure of F{s,p) and P{s,p) as functions of p. 



then F {s,p) is everywhere analytic, except for poles at p = in^ N (with N 
being the natural numbers 1, 2, • • • ), and P {s,p) has poles at p = —in^N. 
In our case k_ = 1 and, as expected, the function F has poles at p = iN 
due to the gamma function^ F (1 + ip). On the other hand, for the case of 
the specific potential (|3H) . one has k+ ^ 0, with the exponential behaviour 
replaced with a power law This is a direct consequence of the fact that 
the metric distance from a point in region Ij„ to the horizon at s — > oo is 
finite, whereas it is infinite going to past infinity at s — > — oo. Therefore, 
the poles at — i^+N become infinitely close and are effectively replaced by a 
branch cut along the negative imaginary p axis, as one can check from the 
analytic expression for P. The features just described, which are shown in 
Figure ^2 are generic and do not depend on the details of the potential, 
but only on the asymptotic behaviour (|33() . In particular, one has the same 
behaviour for fiuctuations of any field. 

Let us now consider a general solution A{s,p) of Since both 

F{s,±p) and P{s,±p) are bases of the solutions of (|32|l. one must have 
that 

Ais,p) = Vf{p)F{s,p) + Uf{p)F{s,-p) 

= Vp ip) P is,p) + Up (p) P (s, -p) . (36) 

The coefficients Vp^p and Up^p can be easily computed by recalling that, 

^The function Jip (—pt) is an entire function of p, as can be seen from the explicit 

power series representation of the Bessel function (z). 
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Figure 12: Figure (a) shows the various coordinate systems in region 
Figure (b) shows the various ways to define a solution A. One either gives 
the functions Vp, Up or the functions Vp, Up, or, alternatively, the functions 
A(so,x), A(so,x) at some fixed time sq along the dotted line. 

given two solutions f,goi (|32)) . the Wronskian 

[/, 9] = fa- 9f 

is independent of s, where the dot denotes differentiation with respect to s. 
In particular, one can easily show that 

[F{s,p),F{s,-p)] = [Pis,p),Pis,-p)] = 2ip. 

Therefore we readily conclude that 

Vp{p) = ^[A{s,p),F{s,-p)] , (37) 

Zip 

Up{p) = -^[Ais,p),F{s,p)] , 
zip 

together with similar equations for Vp and Up. 
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5.4 Behaviour at the Horizons 



The meaning of the coefficients Vp and Up is quite clearly understood by 
considering the value of the field A on the horizons. Since we are in the limit 
s ^ cx), we use the asymptotic form of F (s,p) to conclude that (see Figure 

A (v, u) ^ Vf (v) + Up (u) , (38) 

where 

Vf (v) = J dpVF (p) e-'P" , Uf (n) = j dpUp (p) e^^" . 

To compute the stress tensor due to the fluctuation in A, whose divergence 
usually signals the instability of the Cauchy horizon, we also must consider 
the derivatives of the field with respect to the coordinates (the coordi- 
nates u, V tend to oo on the horizons). Recalling that 

d+ = -e^d^ , d- = -e'^du , 

we conclude that 

d+A{v,u) -e^d^VFiv), (39) 

d-A{v,u) -e^'duUFiu). 

Note that both expressions (|38|) and (|39|) are valid only in the limit (u — > oo, 
V fixed) or (t; ^ oo, M fixed), and are not approximate expressions for large 
but finite u, v. This is because, although the asymptotic form (|34j) is exactly 
valid in the limit s ^ oo, it becomes accurate for s ^ s(p), where the 
threshold value s(p) depends explicitly on p. Given the explicit form of 
F, it is clear that s(p) ~ Inp since only in this case is the argument of 
the Bessel function —pt <S 1. Therefore, given a function A which has 
support at all momenta, one is not uniformly in the asymptotic region for 
all momenta at any large but finite s. Only in the limit s ^ oo we can use 
uniform convergence to deduce the expressions (|nH|) and (jHH]) . in the sense 
just described. 

Finally, going back to the computation of the derivatives 9±A, one re- 
calls from j37j that the problematic limits circ given by lim^ 9+ A and 

limu^oo 5_A. Concentrating on the first one for simplicity, we see that we 
need to consider the expression 

lim 9+ A = lim e" [ dpVp {p) ip e'''^'" (40) 

11— »oo t)— »oo J 

= \ lim e" I dp[A{s,p),F{s,-p)] e-'P\ 

2 f— ^oo J 
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The physical problem we need to address is the following. Let us assume 
that, at some time sq ^ in the past, much before the scattering potential 
V (s) becomes relevant, our field A is in some specific configuration given by 
initial conditions on the value of A and of its time derivative A 

A{so,x) , A(so,x) . 

We will assume that the functions above are localized as a function of x. 
More precisely, we shall demand that A(so,a;) and A{sq,x) are such that 
their evolution would not lead to an inconsistent behaviour on the horizon 
in the Milne wedge of flat space, where V = 0. It is important to realize 
that those are the general physical perturbation, since we are only demand- 
ing that they would be regular in flat space. Let us be more explicit. In 
the absence of the potential V, the solutions of the massless Klein-Gordon 
equations are just A = a{v) + /?(n), and the initial conditions A and A at sq 
are equivalent to giving the functions a and /? of the light -cone coordinates. 
Recall though that the coordinates which are regular across the horizon are 
z^, and therefore we must demand, in order to have a regular perturbation 
in flat space, that a and /? be well-behaved as functions of z^. To see what 
this means in practice, consider a function of z^ which has a nice power 
series expansion oq + aiz^ + ■ ■ ■ around z+ = 0. As a function of v this 
reads a{v) = oq — aie~^ + • • • , and this shows that, aside from the constant 
part, the function a{v) must decay, for v — > oo, at least as fast as e"^. This 
fact imposes constraints on the Fourier transforms of a and /? and, in turn, 
on the functions A(so,x) and A{so,x). In practice, it is sufficient to require 
that the Fourier transforms A(so,p) and A{sq,p) do not have poles in the 
strip |Imp| < 1. Then the problematic limit 



— lim e" I dp 



A (so,p) F {so, -p) - A {so,p) F (sq, -p) 



-ipv 



can be computed by deforming the contour in the lower half of the complex 
plane, and is determined by the poles of F{s,—p). In fact, the leading 
behaviour of the integral is controlled by the first pole at —i, which gives 
an asymptotic behaviour for large v ^ oo of . The limit in the above 
expression then tends to a finite result. 

In order to better understand the above result, let us review the standard 
argument for the instability of the Cauchy horizon, which was followed in 
j54j . The problem is usually posed by giving initial conditions at past null 
infinity (see also Figure ^Ja), by giving the functions 

Vp (v) = J dpVp ip) e-^'P'' , Up {u) = J dp Up [p) e'P"" , 
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on which we impose a regularity constraint (say, for simphcity, a locahzation 
condition hke the one discussed above, now in the coordinates v and u). To 
determine Vp, Up in terms of Vp, Up one uses the expression ()36p for A in 
the equation (|H7|) and obtains that 



Vp{p) \^( A{p) Bi-p) \ ( Vp{p) 
Up{p) ) \B{p) A{-p) ) \ Up{p) 



where^ 



A{p) = ^[P{s,p),F{s,-p)] , 

B{p) = --L[p(s,p),F(s,p)] . 
lip 

The above expression shows that both A {p) and B (p) have a cut along the 
negative imaginary p axis (the limit of the poles at —in^n for — > 0). 
The full analytic structure of A{p) and B{p) is shown in Figure [T3l Since 

ip) = Vp (p) A (p) + Up (p) B {—p) we conclude that the leading pole 
of Vp{p)A{p) at —iK^ — > —iO determines the asymptotic behaviour of the 
integral in H40|) to be e~'^+'". Therefore, the limit limt,^oo diverges, 
thus naively signaling an instability of the horizon. 

The above reasoning depends crucially on the assumption that the func- 
tions Up, Vp have a nice Fourier transform (analytic at least in the region 
|Imp| < 1), so that the perturbation is localized on the null directions at 
past null infinity. If this is the case, though, expression ()36|) shows that the 
function A (s, x) is already delocalized (the Fourier transform has a branch 
cut due to P {s,p)) at any finite time s <C —1 much before the scattering po- 
tential, and therefore will concentrate on the horizons and create an infinite 
stress tensor. This is clearly not the type of perturbation we want to focus 
on, which should be localized in space before they hit the potential V. Note 



**The coefficients A (p) and B (p) satisfy, generically, the conjugation relations A (p) — 
A {—p) and B (p) — B (— p), and the unitarity relation A (p) A (— p) — B (p) B (— p) — 1. 



Using 13511 . we have the explicit expression 



sinh (vrp) F (1 + ip) ' 



2 V2/ sinh(7rp) r(l - ip) 

Moreover, in general, the function A{p) will have zeros on the positive imaginary p-axis 
whenever the potential V has a bound state. The explicit expressions above show that, 
for our particular potential, there are no bound states. 



34 



A(p) 



3i « 
2i it 
i n 



B(p) 




-1 



-2i 



-3i 



Figure 13: The analytic structure of A{p) and B{p). 



that, at any finite time s, we cannot use, in equation the asymptotic 
form of P{s,p) for all values of p. In fact, the explicit form of the the function 
P shows that one is in the asymptotic regime for \p\ ^ s~^. Therefore, the 
crucial low momentum modes which determine the asymptotic behaviour of 
the wave function are never in the asymptotic region for any finite time s. 
We then consider the requirement of localization of A (s, x) at times s <C — 1 
to be the correct and physically relevant boundary condition to study issues 
of stability of the geometry. The reasonings of this section then show that, 
in this sense, the Cauchy horizon in perfectly stable to small perturbations. 

5.5 Discussion 

Let us conclude this section with some comments on future research. 

We have shown that, in the orientifold cosmology, one can define particle- 
like perturbations which do not destabilize the Cauchy horizon, and there- 
fore that the solution is stable against small variations of the background 
fields. All the work is done at the classical linearized level. The most press- 
ing question for issues of stability is the study of the quantum stress-energy 
tensor for the field ^. This is a non-trivial problem due to the non-minimal 
coupling to the metric and the dilaton, and only some results are known 
in the literature (see |5n[ I51j and references therein). On the other hand, 
the problem might be tractable since the exact wave-functions are known 
in this case. 
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6 Conclusion 



In this paper wc investigate the classical stability of a two-dimensional ori- 
entifold cosmology, related to a time-dependent orbifold of fiat space-time. 
We saw, with a specific counter example, that the instability argument of 
Horowitz and Polchinski is not valid once the time-like orbifold singularity 
is interpreted as the boundary of space-time. In this specific example, we 
consider an exact shock-wave solution of two-dimensional dilaton gravity 
that uplifts to a distribution of matter in the three-dimensional covering 
space. According to Horowitz and Polchinski this distribution should in- 
teract gravitationally creating large black holes. Indeed, there are closed 
time-like curves in the covering space geometry of the shock-wave, signal- 
ing the three-dimensional gravitational instability. However, such CTC's 
arc not present if we interpret the singularity as a boundary of space-time, 
and accordingly excise from the geometry the region behind it. 

The other stability problem addressed in this work is related to cosmic 
censorship. The presence of naked singularities with a Cauchy Horizon has 
lead many people to believe that such singularities never form. Indeed, as 
for the Reissner-Nordstrom black hole one expects the Cauchy horizon to 
be unstable when crossed by matter. We analyze the propagation of a scalar 
field coupled to gravity in the geometry, and show that any localized fluc- 
tuation at some finite time in the far past will not destabilize the horizon. 
The existence of the time-like singularity does not imply the break-down 
of predictability because of the conjectured duality between the singularity 
and orientifolds of string theory. To leading approximation, the effect of the 
orientifolds is to enforce a boundary condition on the fields, which deter- 
mines uniquely their evolution. A more complete quantum understanding 
of this system is desirable. 
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